IN THE CLAIMS 



1 . (Currently amended) A method for the accurate determination of van der Waals 
parameters for high-precision determination of crystal structures and/or energies, comprising the 
steps of: 

[[-]]numerically simulating at least one crystal structure based on density functional 
theory (DFT) calculations combined with a potential energy term representing Van der 
Waals interactions; 

[[-]]providing reference data containing accurate information about said at least one 
crystal structure; 

[[-]] defining a deviation function (F) quantifying a deviation between said reference data 
and said at least one simulated crystal structure; 

[[-]] fitting at least one parameter of said van der Waals potential term in such a way as to 
minimize said deviation function (F); and 

[[-]]obtaining the accurate van der Waals parameters from the best fit. 

2. (Original) A method according to claim 1 , characterized in that said van der Waals 
potential term is defined as: 



wherein fA,B (rA,B) is a damping function and the sum runs over all pairs of interacting atoms, and 
that said fitting step comprises fitting said damping function. 

3. (Original) A method according to claim 2, characterized in that said damping 
function is defined as 
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and that said fitting step comprises fitting the parameter r a,b, and/or the parameter n and/or the 
parameter c. 

4. (Currently amended) A method according to claim 2-&¥-^, characterized in that said 
fitting step furthermore comprises fitting said coefficient C6,a,b- 

5. (Currently amended) A method according to any of the preceding claims claim 1 . 
characterized in that said reference data are theoretical data obtained by Hartree-Fock 
calculations or QuMitum Monte Carlo simulations. 

6. (Currently amended) A method according to any of the preceding claims claim L 
characterized in that said reference data are experimental low-temperature crystal structure data. 

7. (Original) A method according to claim 6, characterized in that said crystal structure 
data are obtained by X-Ray or neutron scattering. 

8. (Currently amended) A method according to claim 1 for the accurate determination of 
crystal structures and/or energies, comprising the steps of: 

[[-]]providing a rough estimate model of at least one crystal structure; 

[[-]]numerically simulating said at least one crystal structure based on density fimctional 
theory (DFT) calculations combined with a potential energy term representing Van der Waals 
interactions; and 

[[-]] obtaining said at least one crystal structure and/or its energy as a result of said 
numerical simulationJX,]] 

characteriz e d in that said van dor Waals pot e ntial term is obtained by the method according to 
any of claims 1 to 7 . 

9. (Original) A method according to claim 8, characterized in that a plurality of 
polymorphic crystal structures are determined and ranked according to their respective energies. 
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1 0. (Original) A method for the efficient numerical optimization of a molecular crystal 
structure using an advantageous crystal coordinate system, comprising the steps of: 

[[-]]providing a starting crystal lattice described by an initial coordinate system 
comprising lattice parameters and atomic positions in said crystal; 

[[-]]defining a so-called natural coordinate system and representing said starting crystal 
lattice in said natural coordinate system, said natural coordinate system comprising: 

[[•]]first coordinates describing symmetry-allowed lattice changes and defined in 
such a way that changes of said first coordinates do not cause changes of the 
molecular geometry or a rotation of molecules with respect to each other and 
leave fractional coordinates of molecular centres constant; 
[[•]]second coordinates describing symmetry-allowed translations of said 

molecules in said crystal; 
[[•]]third coordinates describing symmetry-allowed rotations of said molecules in 
said crystal; 

[[•]]fo\irth coordinates describing symmetry-allowed changes of the molecular 
geometry; 

[[-]] transforming coordinates from said natural coordinate system to said initial 

coordinate system; 

[[-]]calculating the lattice energy and energy derivatives with respect to said initial 
coordinate system; and 

[[-]]transforming said energy derivatives fi-om said initial coordinate system to said 
natural coordinate system, 

wherein a minimization algorithm is used for minimizing said lattice energy with respect to said 
natural coordinate system. 

11. (Original) A method for the energy ranking of polymorphic crystal structures, 
comprising the steps of: 

[[-]]providing rough estimate models of each of said crystal structures; 

[[-]]numerically simulating each of said crystal structures based on density functional 
theory (DFT) calculations combined with a potential energy term representing Van der Waals 
interactions 
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[[-]]obtaining accurate crystal structures and energies as a result of said numerical 
simulation; and 

[[-]]ranking said accurate crystal structures according to their respective accurate 
energies. 

12. (Original) The method according to claim 1 1 , characterized in that the crystals are 
crystals of pharmaceutical compounds. 

1 3 . (Original) The method according to claim 1 2, characterized in that it is applied to 
identify the most stable polymorphic form of a pharmaceutical compound. 

14. (Currently amended) The method according to any of claims 1 1 to 1 3 

7 characteriz e d in that it us e s a method for the effici e nt optimization of the molecular crystal 
structur e according to claun 10 . 
claim 1 1 including: 

providing a starting crystal lattice described by an initial coordinate system comprising 
lattice parameters and atomic positions in said crystal: 

defining a so-called natural coordinate system and representing said starting crystal 
lattice in said natural coordinate system, said natural coordinate system comprising: 

first coordinates describing symmetry-allowed lattice changes and defined 
in such a way that changes of said first coordinates do not cause changes of 
the molecular geometry or a rotation of molecules with respect to each other 
and leave fractional coordinates of molecular centres constant: 

second coordinates describing symmetry-allowed translations of said 
molecules in said crystal: 

third coordinates describing symmetry-allowed rotations of said molecules 
in said crystal: 

fourth coordinates describing symmetry-allowed changes of the molecular 

geometry: 

transforming coordinates from said natural coordinate system to said initial coordinate 

system: 
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calculating the lattice energy and energy derivatives with respect to said initial coordinate 
system: and 

transforming said energy derivatives from said initial coordinate system to said natural 
coordinate system. 

wherein a minimization algorithm is used for minimizing said lattice energy with respect 
to said natural coordinate system. 

15. (Currently Amended) A comput e r program product comprising computer readable code 
for e nabling a comput e r to perform a m e thod according to any of the pr e c e ding claims when said 

code is ex e cuted. The method according to claim 1 wherein the steps for determining the van der 
Waals parameters are performed by a computer program comprising computer readable code 
executable by a computer. 

1 6. (New) The method according to claim 8 including providing a computer program 
comprising computer readable code executable by a computer that determines the crystal 
structures and/or energies. 

1 7. (New) The method according to claim 1 0 including providing a computer program 
comprising computer readable code executable by a computer that numerically optimizes the 
molecular crystal structvire. 

1 8. (New) The method according to claim 1 1 including providing a computer program 
comprising computer readable code executable by a computer that energy ranks the polymorphic 
crystal structures. 
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